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I.  INTRODUCTION 


Many  laser  systems  operating  in  the  atmosphere  can  be  degraded 
atmospheric  turbulence.  For  example,  target-illumination  systems,  communication 
systems,  radar  systems,  and  others  may  be  severely  affected  by  turbulence  effects 
including  beam  scintillation.  For  this  reason  a  considerable  effort  has  been  made 
to  understand  and  correct  for  these  turbulence  effects^^^.  These  atmospheric 
correction  efforts  are  based  on  determining  the  phase  distortions  introduced  by 
turbulence  and  correcting  them  at  the  laser  ap«*ture.  The  theory  that  has  been 
developed  along  these  lines  to  date  has  been  in  the  regime  of  weak  turbulence  only, 
where  the  beam  wanders  as  a  whole  and  does  not  break  up  into  multiple  patdies  or 
blobs.  Under  this  condition  a  correction  for  wavefront  tQt  leads  to  a  considerable 
improvement  of  on-axis  irradiance. 

In  this  paper,  we  are  concerned  with  the  effectiveness  of  this  tilt  correction 
as  the  effects  of  turbulence  become  more  pronounced  (i.e.,  where  in  fact  the  beam 
breaks  up  into  multiple  patches).  To  this  end  we  develop  formulae  for  the  mean 
square  beam  wander  angle  of  a  focused  laser  beam  transmitted  through  a 
homogeneous-isotropic  (e.g.,  constant  altitude)  atmosphere  imder  both  weak  and 
strong  turbulence  conditions.  We  find  that  tat  correction  is  of  marginal  benefit  (at 
best)  under  strong  turbulence  conditions.  Since  atmospheric  phase  aberrations 
must  be  corrected  at  the  laser  aperture,  the  phase  of  a  spherical  wave  emitted 
from  a  source  in  the  focal  plane  of  the  laser  is  measured  and  the  tilt  of  this  wave  is 
determined.  The  "negative"  of  this  tat  is  then  used  as  the  correction.  It  is  shown 
that  differences  exist  between  the  wander  angle  of  the  transmitted  laser  beam  and 
the  tat  of  the  received  spherical  wave.  These  differences  exist  since  it  is  not 
possible  to  collect  the  entire  field  due  to  the  point  source  in  the  laser  aperture 
plane.  However,  in  practice  these  differences  appear  to  be  small  indicating  that 
the  method  of  tilt  correction  is  useful. 

In  this  report  we  will  use  the  centroid  as  a  definition  of  wander  angle  or  tilt. 
Since  most  adaptive  optical  systems  measure  "phase"  and  define  the  tat  in  terms  of 
a  least  squares  fit  over  the  aperture,  it  is  of  interest  to  determine  the  correlation 


of  the  tilt  defined  by  the  two  methods.  We  have  only  been  aUe  to  do  this  in  the 
weak  turbulence  regime  where  the  correlation  between  the  two  methods  is  found  to 
be  nearly  100%. 

Engineering  formulas  are  also  developed  for  quick  calculation  of  wander 
angle  and  peak  on-axis  irradiance  of  the  transmitted  beam  for  which  a  tat 
correction  has  been  made.  These  formulas  are  simple  enough  that  hand  calculation 
suffices. 

In  the  next  section  we  will  define  the  regimes  of  weak  and  strong  turbulence 
and  briefly  describe  the  results  of  this  work  for  those  who  are  not  interested  in  the 
theoretical  details.  In  Section  III  we  will  derive  general  expressions  for  the  2-axis 
mean  square  wanda*  angle  for  both  the  transmitted  laser  beam  and  the  received 
spherical  wave.  Specific  derivations  for  the  two  cases  will  then  foUow  in  Sections 
IV  and  V.  In  Section  VI  we  will  derive  the  correlation  function  for  tilt  defined  in 
terms  of  the  centroid  and  the  tilt  defined  in  the  least  squares  sense.  Finally  in 
Section  VII  we  will  present  engineering  formulas  and  comment  on  their  applica¬ 
bility. 
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n.  DEFINITIONS  AND  RESULTS 


An  optical  beam  interacts  with  atmospheric  refractive  inhomogeneities  of 
all  scales.  Turbulent  eddies  having  characteristic  dimensions  greater  than  the 
beam  diameter  lead  to  refractive  effects  (beam  wander);  these  large-scale 
inhomogeneities  lead  to  net  wavefront  tilt.  Conversely,  eddies  having  character¬ 
istic  dimensions  less  than  the  beam  diameter  lead  to  diffractive  effects  about  the 
instantaneous  center  of  energy.  The  small  eddies  cause  the  primary  amplitude 
effects  on  the  wavefront.  These  effects  become  more  pronounced  as  the  strength 
of  turbulence  and/or  propagation  range  increases.  A  quantitative  measure  which 
separates  the  regions  of  weak  (amplitude  effects  negligible)  from  strong  turbulence 

(strong  amplitude  effects)  is  given  by  the  first  order  Rytov  log-amplitude  variance 

(2) 

of  a  spherical  wave  which  for  constant  altitude  propagation  is  given  by 


=  0.124 


(1) 


2 

where  k  =  2ir/X,  X  is  the  optical  wavelength,  z  is  the  propagation  range,  and  is 

the  index  of  refraction  structure  function.  As  is  well  known,  the  log-amplitude 

variance  saturates  to  a  value  close  to  unity  as  the  propagation  range  and/or 

increases.  Equation  (1)  does  not  behave  in  this  manner;  however,  the  propagation 

2 

range  dividing  weak  from  strong  turbulence  can  be  defined  by  setting  =  1.  This 
range  is  referred  to  here  as  the  saturation  distance  and  is  given  by 


z 

s 


0.124 

n 


6/11 


(c2)6rtV/ll 


(2) 


2 

A  plot  of  this  saturation  range  vs  wavelength  with  as  a  parameter  is 
shown  in  Fig.  1.  As  is  seen  from  this  figure,  amplitude  effects  are  much  more 
pronounced  for  short  wavelength  lasers.  A  measure  of  the  phase  degradation  of  a 
wave  is  given  by  the  long  term  lateral  coherence  length  of  a  spherical  wave, 
defined  by^^^  the  detector  separation  at  which  the  mutual  coherence  function  M(p^) 
—  €  ,  i.c. 
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(3) 


2  2 

p  (z)  =  (0.545  V.  C  z) 
n 

2 

We  have  plotted  vs  wavelength  with  as  a  parameter  in  Fig.  2.  As  can  be 

seen  from  this  figure,  P^^^  is  considerably  smaller  than  the  aperture  diameter  for 
most  cases  of  interest.  We  will  also  define  two  dimensionless  parameters  in  which 
the  results  of  this  paper  are  e3?)ressed.  These  are  the  ratio  of  the  laser  or 
receiving  aperture  diameter  to  and  the  normalized  propagation  range,  i.e., 

c  =  D/p  (z  );  z  =  z/z  (4) 

OS  n  s 

2 

Finally,  to  complete  our  definitions,  we  present  a  plot  of  vs  altitude  in  Fig.  3 
which  is  useful  in  calculating  various  cases  for  the  wander  angle.  The  model  shown 


=  5  X  1  <  H  <  T-2000 

n 

T<  H  <  « 


=  8  X  10"^^  T-2000  <  H  <  T  (5) 

n 


2  -2/3 

where  the  units  of  H  are  in  meters,  C  is  m  ,  and  T  is  the  altitude  of  the 

2  ^  (5) 

tropopause.  Other  models  for  exist,  notably  those  developed  by  Hufnagel. 

With  these  definitions  we  now  present  the  results  of  the  present  study.  The 

2-axis  mean  square  wander  angle  of  a  focused  laser  beam  under  weak  and  strong 

turbulence  conditions  vs  the  aperture  otascuration  ratio  b  is  (the  1-axis  value  is 

equal  to  one-half  the  corresponding  2-axis  value) 


^''T'^weak 


400  €®^^z„ 

«  /  x^'''M,(x)e 


h  2/3  4B^(Dx) 

f  x2'^M.(x)e  dx 


(6) 
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Fig.  2.  Long  Term  Lateral  Coh 


400  1 

^^T>strong  =  - ^ - 72  I  [l  -  (x)expr-2(€x)®/\  Jdx 

(7) 

where  D  is  the  aperture  diameter,  M  .(x)  is  the  aperture  mutual  coherence  function 

(6)  ® 
given  by'  ' 


H(x)  =  1  x>0 


=  0  x<  0 


and  is  the  log-amplitude  correlation  function  of  a  spherical  wave.  Under  all 
turbulence  conditions  Clifford  has  shown^^^ 


I 


If  we  apply  asymptotic  limits  to  both  the  weak  and  strong  turbulence  conditions  we 
can  obtain  a  connecting  formula  which  appears  to  be  valid  for  all  turbulence 
regimes.  The  resulting  expression  for  the  2-axis  mean  square  transmitted  wander 
angle  is 


<ep  = 


400  z^e 


5/3 


AW) 


(kD)' 


11/3 


1  + 


0.5z^^^^^l-(J^) 


(10) 


where 


A(i)  =  /  x^^^M^(x)dx 


(11) 


=  0.275  -  0.989  6 


2.29 


0  <  d  <  0.5 


For  the  received  spherical  wave  case,  the  results  are  essentially  identical  to 
the  transmitted  laser  beam  case,  except  under  strong  turbulence  conditions  (z^z^) 
for  which 


<0p  =  G  (12) 

where 


00 

1 

F 

\  F/ 

(2i+l)!!  ’ 

L 

i=0  -1 

48/16/5) 

■  / 

5)r 

and  T^is  the  gamma  function.  Note  that  F  is  the  normalized  variance  of  received 
power  collected  by  the  aperture  due  to  the  spherical  wave.  This  variance  is 
generally  very  small  for  strong  turbulence  (due  to  aperture  averaging)  and  an 
asymptotic  expression  exists  for  G: 


=  1  +  3F  +  15F^  + 

asym 


(13) 


Noting  that  the  first  order  correction  to  G  decreases  as  z 
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and  that  similar 
2> 


contributions  were  neglected  in  deriving  the  strong  turbulence  limits  of  we 

2  2  * 
can  conclude  that  to  a  numerical  factor  of  the  order  of  unity. 

As  indicated  above,  for  adaptive  optical  systems  there  is  interest  in 

determining  the  correlation  between  the  wander  angle  defined  by  the  centroid  ^ 

and  the  tilt  defined  by  a  least  squares  fit  of  the  phase  over  the  aperture,  as 
P  /Q\  2 

defined  by  Fried  .  Note  ^  is  used  to  obtain  For  weak  turbulence 

conditions,  i.e.,  z  <  we  find  that  the  correlation  is  given  by 


e  e 

p  c 


(i)  = 


■  ^ 


I 


x®'^^M^(x)dx 


N 


1/2 


ryx2«M,(x)dxj‘^2 


(14) 


where 


and 


= 


Hr, 

X  /74/3)  180 


r(ll/6)  47311 


-i/ 


l*d 

2 


1  +<» 


23/3J . 


11/3 


Id. 

2 


rT(^,  5,  a,  ^)d^ 


(15) 


T(^  e,  a,  ?■)  =  d*a  +1 


j  0  sino  -  0  +  j  ^  sin  ? 


+  16  sin^a  /cot  ?  -  cot  «  +  cqtV-.cot^a^ 


64  4  .  3 
-  ^  sin  a 


(csc^?  -  csc^a)  +  sin  a(cot  0  -  cot  a) 
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with  ?  =  a  -  and  «,  0  being  defined  in  Eq.  8.  For  no  obscuration 

Cg  n  (0)  =  0.9919  while  at  d  =  0.7,  C  «  ^  (0.7)  =  0.9602,  indicating  a  high  degree 

*c  p 

of  correlation  for  weak  turbulence  conditions  for  all  cases  of  practical  concern. 
Due  to  mathematical  complexity  we  have  not  been  able  to  determine  the  results  in 
the  strong  turbulence  limit.  However  we  expect  the  correlation  to  decrease  as 
turbulence  increases  since  the  phase  becomes  uniformly  distributed  over  2ir.  A 
measurement  of  phase  using  an  interferometer  may  therefore  be  meaningless  under 
strong  turbulence  conditions. 

Finally,  an  engineering  formula  for  is  given.  This  quantity  is  the  ratio 
of  the  peak  corresponding  on-axis  irradiance  with  tat  correction  to  the  correspond¬ 
ing  on-axis  irradiance  in  the  absence  of  turbulence.  We  find  that 


where 


'rel 


(16) 


«vac  =  2X/xD, 


=  f? 

^vac 


i(<j)  = 


400z  e 
n 


5/3 


9ir(l-i^)(kD)^ 


-*>B  .  z 

n 


(17) 


«g  =  2.51  X  10"^  , 
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T  is  the  atmospheric  attenuation  coefficient  due  to  all  inherent  linear  extinction 

characteristics  of  molecules,  aerosols,  clouds,  etc.,  is  the  2-sigma  mean  square 

jitter  an^e,  m  is  the  aperture  diffraction  limit  performance  factor  (m  is  unity  for 

a  perfect  aperture,  it  is  usually  bigger  than  1  due  to  imperfections  in  the  aperture 

2 

or  aperture  illumination),  and  the  laser  beam  mean  square  wander  angle,  is 
defined  in  Eq.  10.  In  the  next  section,  the  mean  square  angle  of  arrival  will  be 
defined  in  terms  of  the  centroid  and  the  general  expression  derived  for  both  the 
transmitted  laser  beam  and  received  spherical  wave  cases. 
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m.  FORMULATION  OF  WANDER  ANGLE 


We  are  interested  in  treating  two  wander  angle  problems.  The  first  is  the 

determination  of  the  mean  square  wander  angle  of  a  transmitted  laser  beam  in  the 

far  field  or  focal  region  of  the  laser  aperture  (Fraunhofer  regime).  The  second 

problem  is  the  determination  of  the  mean  square  angle  of  arrival  of  a  spherical 

wave  generated  by  a  point  source  in  the  laser  target  plane  and  observed  by  the 

same  aperture  generating  the  laser  signal.  These  problems  are  treated  statistically 

with  anglular  fluctuations  caused  by  the  random  changes  in  the  index  of  refraction 

of  the  atmosphere  due  to  temperature  variations.  It  is  assumed  that  the  index  of 

refraction  fluctuations  are  isotropic  and  homogeneous  and  can  be  described  by  the 

Kolmogorov  spectrum.  The  strength  of  turbulence  is  described  by  the  index 
2 

structure  constant  which  is  assumed  to  be  independent  of  propagation  distance. 

The  extension  to  non-uniform  turbulence  conditions  is  straightforward. 

For  both  cases  under  consideration  the  optical  field  in  angular  coordinates 

(2) 

(a,  0)  in  the  focal  plane  can  be  represented  as  ' 


U(«,  U)  =  271  ^-R-^ y ^  ^^^dr/d^  (18) 

where  U^(J7,  |)  is  the  field  at  the  position  (H,  0  the  aperture  plane  due  to  a 
spherical  wave  from  a  point  source  located  in  the  laser  target  plane  (note  that  the 
point  source  is  located  at  a  position  denoted  by  (a,  p)  in  the  transmission  case), 
W(J7,  ^ )  is  the  aperture  weighting  function,  k  =  2ir/X,  R^  is  the  far  field  "focal" 
distance,  and  S  represents  integration  over  the  clear  aperture  area.  In  the 
transmission  problem  the  weighting  function  includes  the  field  distribution  of  the 
laser  beam  incident  on  the  aperture.  For  both  cases  under  consideration  we  will 
assume  that  the  weighting  function,  except  for  obscuration,  is  uniform,  i.e.,  W(n,4  ) 
=  1.  Non-uniform  weighting  introduces  extra  terms  in  the  analysis  which  will  be 
presented  but  not  carried  as  we  complete  the  analysis.  The  angular  distribution  of 
intensity  is  given  by  the  product  U(a,  P)U '  (a,  Ph 
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l(a,  0)  =  ur  =  ^')W(r7,Ow''(n',  O 

g-il<[a(J7-;7)  +  ^gj 

The  angular  coordinates  of  the  center  of  gravity  (centroid)  in  the  focal  plane 
for  both  cases  of  interest  is  defined  as 


JJ"  I(a,0)a6a60  JJ"  l(a,P)0dadP 
-00  „  -00 

/r"  ’ 

JJHa,0)dad0  JJl{a,0)dad0 


where  the  integrations  are  carried  out  over  the  entire  focal  plane.  Let  us  first 
evaluate  the  integral  in  the  denominator. 

For  the  transmission  case  the  denominator,  by  energy  considerations,  is  the 
total  power  emitted  by  the  aperture.  For  the  receiver  case 


0)dadP  = 


fffk,. 


4ir2R2  -g-g 


dad^ 


00 

-  00 

=  |UQ(^,^)|^|w(77,4)l^d7?d| 

E  (21) 

This  is  the  power  intercepted  from  the  point  source  by  the  aperture.  This  quantity 
is  not  a  constant  but  varies  in  a  statistical  manner  and  hence  must  be  included  in 
any  ensemble  averaging  process.  The  evaluation  ol  the  numerator  follows  in  a 
similar  manner: 
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/“““  ,0)da 


JJJJ 

”  EE 


jOOfjTwtt, 


VVi"(rj  ,Odf7d^d;7d^ 


wjf. 


For  the  receiver  case  U  and  U  '  are  independent  of  (a,P)  and  this  term  can  be 

o  o 

removed  from  the  integration.  For  the  transmission  case  an  argument  can  be  made 


that  if  k^a^+p  D  ^  1  then  the  exponential  terms  oscillate  very  rapidly  giving 
zero  contribution.  This  means  that  significant  contributions  to  the  integral  occur 

.  -  - >4^  +  P^<  ^ 


only  for  angular  coordinates  such  that 


However, 


varies  on  scales  of  the  order  of  -  where  p.  is  the  lateral 
o  o  /,»  o 

coherence  length.'  For  D  >  p^,  which  is  the  case  of  interest  here, 
ctm  be  considered  independent  of  a,  P.  This  again  implies  that 
these  terms  can  be  held  constant  for  integration  over  a,  P;  i.e. 


JJ aI(a,p)dadP  =  U^(rj,4)W(rj,() 


dlv*(n,Owh,()] 


drjd^  (23) 


This  is  obtained  by  noting  that' 


. 

y-iua  . 
ae  da  - 


2ird'(u)i 


At  this  point  we  drop  the  aperture  function  W(>7,0  assuming  it  has  a  zero  gradient 
and  noting  that  its  inclusion  results  in  an  extra  term  in  the  numerator. 

The  mean  square  wander  angle  for  the  transmission  case  is  then  equal  to 


<»P  -  <‘■>1 


1 

^  -  -  3 
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where  P  is  the  power  transmitted  by  the  aperture,  and  <>  indicates  the 
ensemble  average.  The  corresponding  mean  square  wander  angle  in  the  receiving 
case  is  equal  to 


<«R> 


Uo(r)[^Jr)][£Ujj(r')]uJr') 

(yu^(2)U^(2)dp)^ 


^  drdr ' 


We  consider  these  cases  separately  below. 


-20- 


n 


IV.  FAR-FIELD  MEAN  SQUARE  WANDER  ANGLE  OF  A  TRANSMITTED 

LASER  BEAM 

To  obtain  the  mean  square  wander  angle  for  the  transmitted  U  er  beam  the 
ensemble  average  of  U^(r)[7U^(r)]  •  (7U  (r»Uj^(r')  is  required.  This  average  can  not 
be  determined  over  all  turbulence  regimes  since  the  field  statistics  are  not  known 
for  the  intermediate  turbulence  strength  regime  (i.e.,  ~  1).  Instead  we  consider 

the  regime  of  weak  turbulence  conditions  (z^  «  1)  for  which  the  field  statistics 
are  known  to  be  log-normal  and  the  case  of  very  strong  turbulence  conditions  (z^ 
»  1)  for  which  the  field  statistics  is  known  to  be  normal.  We  then  use  a  smooth 
connecting  formula  to  obtain  results  which  we  postulate  are  valid  over  all 
turbulence  conditions  regimes. 


A.  WEAK  TURBULENCE  RESULTS 


In  the  weak  turbulence  conditions  regime  the  field  statistics  are  log- 
(2) 

normal.  Therefore  we  let 


X(r)  +  iS(r) 
e  —  — 


(26) 


where  X(r)  is  the  log  amplitude  and  S  the  phase  of  the  spherical  wave  from  a  point 
source  located  in  the  target  plane.  The  required  ensemble  average  is  then 
rewritten  as  follows; 


2  3 

13-*^ 

where  we  have  taken  the  gradients  outside  the  ensemble  average  by  the  use  of  the 

limit  operation.  The  validity  of  this  approach  will  be  demonstrated  in  the  section 

on  strong  turbulence.  Writing  the  field  in  terms  of  the  log-amplitude  and  pheise, 

the  fourth  order  correlation  function  in  U^  becomes 

o 
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(28) 


r(ri,r2.r3’Li)  =  <Uo^Ii)Uo^j:2^U  Jr3)U^(r4)> 

=  <exp[X(rj)  ■*’X(r2)  ■'■X(p3)  ■'■X(p4)  +  ®(£^) 

-  iS(r2)  +  'S(r3)  -  iS(r^)]>. 

For  a  XandS  which  are  normally  distributed,  it  can  be  shown  to  second  order  (in 
the  exponent)  in  the  index  fliwtuations  that  <e*^^  =  expKi//^  +  | )^^1. 
Therefore  the  fourth  order  correlation  is  given  by 

r(li.  12, 13,  P4)  =  exp  I  -  i  |^I^(|ri-r2|)  +  D^(|r3-r4l)  +  D^(|rj-r4|)  +  D^(|r2^3|) 

■  2Bj^(|rj  -£3!)  +  2^(|r2-r4l ) 

+  2i«SX(|rj-r3|)>  -  <SX(lr2:r4|)>j  ,  (29) 

where  is  the  wave  structure  function,  is  the  log-amplitude  correlation 
function,  and  SX  is  the  correlation  between  the  phase  and  log-amplitude/^^  Taking 
the  gradients  with  respect  to  coordinates  2  and  3  and  performing  the  limit 
operations  we  obtain  that 

lim 


If  we  substitute  the  resiite  of  Eq.  30  into  Eq.  27  and  then  into  Eq.  24,  and 
'2  '2  2 

neglect  the  terms  [B^]  and  [SXl  in  comparison  to  7  Dj^  (see  Appendix  A)  we 
obtain  the  resulte  for  the  mean  square  wander  angle  of  a  transmitter  laser  beam  as 
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where  Mj(x)  is  the  normalized  overlap  integral  (MTF)  due  to  the  aperture  (Eq.  8) 

and  is  the  uniform  irradiance  at  the  aperture.  Since  all  propagation  paths  are 

far  above  the  ground,  we  neglect  the  effect  of  the  outer  scale  on  the  wave 

structure  function  D^.  (If  the  beam  diameter  becomes  equal  to  or  greater  than  the 

size  of  the  outer  scale  then  the  outer  scale  effects  must  be  included  in  D^.  The 

outer  scale  is  approximately  equal  in  magnitude  to  the  height  above  the  ground  up 

(3) 

to  several  hundred  meters.)  Then 

=  ,  (32) 

where  is  given  in  Eq.  3.  Taking  the  Laplacian  of  j  D^(|p|)  yields 


9p  5/3pl/3 
o 


Using  Eq.  4  in  the  form 


p;5/3(z) 


,5/3^  D-5/3 
n 


yields 


<4>  = 


400e^^^z 


9ir(kD)^(l-d^) 


1 

«  2  0 


2/3 M  (  \ 

X  M.  (x)e  dx, 

0 


where  (x)  and  are  defined  in  Eqs.  8  and  9,  respectively.  For  weak  turbulence 
is  very  small  and  the  exponential  factor  in  Eq.  35  can  be  neglected.  We 
will  include  this  factor  for  more  exact  calculations  but  fca*  the  engineering 
formulas,  developed  below,  this  term  is  not  necessary. 

The  final  result  under  weak  turbulence  conditions  for  the  2-axis  mean  square 
wando*  angle  for  a  transmitted  beam  is 


a.<4>  =  ^ 


(kD)^(l-6^) 


X 

J  x^^^M^(x)dx. 
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B.  STRONG  TURBULENCE  LIMIT 

The  principal  difference  between  the  weak  and  strong  limits  of  turbulence  is 
that  under  strong  turbulence  the  fields  are  statistically  normal  while  in  the  weak 
regime  they  are  statistically  log-normal.  This  leads  to  significant  difference  when 
treating  the  ensemble  average  expressed  in  Eq.  24.  We  will  perform  the  ensemble 
average  in  two  ways.  The  first  follows  the  technique  used  in  Eq.  27  of  using 
limiting  operations  and  taking  the  derivatives  outside  the  ensemble  average  while 
the  second  takes  the  derivatives  of  U^  before  taking  the  ensemble  average.  The 
two  results  are  identical  verifying  the  first  technique  used  here  and  in  the  previous 
section.  Following  the  method  used  in  Eq.  27,  the  determination  of 
£2’ -3’ ^  required  under  strong  turbulence  conditions.  For  strong 
turbulence  conditions  the  field  statistics  are  jointly  normal.  Hence, 

where  uncompensated  phase  terms  have  been  neglected  since  they  wQl  average  to 

(9) 

zero  very  rapidly  due  to  the  high  optical  frequency.  Fante  has  shown  the 

existence  of  a  long  range  correlation  tail  besides  the  terms  given  in  Eq.  37.  We 

2 

will  show  that  this  tail  does  not  contribute  significantly  to  <0rj,>  for  large  in 
appendix  B.  For  strong  turbulence,  it  has  been  shown  that  the  second  order 
correlation  function  is  the  same  as  for  weak  turbulence,  i.e., 

rti.Ij)  =  <U„(r,)uJr2)>  =  expl-^  (38) 

which  implies  that  the  fourth  order  correlation  function  for  strong  turbulence  is 
given  by 
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(39) 


f'^ly  l2*  Is’  £4^  “  2  "  2 

*  «»Pl-  7  -  5  V'i2--E3'’'- 

Performing  the  gradient  and  limiting  op«*ations  yields 

lim  J?2'-^3  -1’ -2’ -3’ “  2  ^ ~  2 
^1— .Z9— .£ 

lU^-r 


(40) 


where  P  =  r  -  r'.  Substituting  this  result  into  Eq.  24,  performing  the  integration 
over  the  sum  coordinate  and  expressing  the  results  in  terms  of  e  and  yields  the 
2-axis  mean  square  wander  angle  of  a  transmitted  laser  beam  in  strong  turbulence 
is  given  by 

„  400€^/^z  } 

B=<0_>  = - 2-^  y  X  '  [l-(fx)^'^z JM  (x)exp[-2(€x)®^^z  ]dx. 

*  oo^n  ^0  " 

9ff(kD)''(l-i^)  “ 

(41) 


This  is  the  main  result  of  this  section. 

Note  that  in  very  strong  turbulence  an  asymptotic  formula  may  be  derived 
from  Eq.  41.  For  this  case,  the  scale  on  x  is  very  small  due  to  the  exponential 
factor;  therefore,  M^(x)  may  be  replaced  by  M^(0)  =  and  the  limit  on  the 

integral  extended  to  infinity: 


^4\sym. 


20 


3(kD)2(l-d2) 


ft. 


y/2)c 


■ydy  =  ^  (kD)"2(l-i‘) 


2v-l 


(42) 


This  impUes  that  in  very  strong  turbulence  the  rms  wander  angle  of  the  transmitted 
laser  beam  tends  to  a  constant  which  is  of  the  same  order  as  the  diffraction  limited 
angular  spot  size  of  the  aperture. 
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We  now  return  to  a  consideration  of  the  second  method  of  determination  of 
the  ensemble  average  expressed  in  Eq.  24.  First  we  note  that  the  field  may  be 
expressed  in  terms  of  its  real  and  imaginary  parts  (X,Y)  which  are  statistically 
independent  and  distributed  normally  with  zero  mean  and  the  same  variance.  The 
ensemble  average  becomes 

e  =  <lX(r)  +  iY(r)][£X(r)-i7Y(r)l  •  [VXCr')  +  i7Y(r')][X(r')-Y(r')]> 

=  (l|£j.[x\)  +  Y^(r)l  -  i[X(r)7pY(r)-Y(r)£pX(r)]|.||7p'[X^(r>Y^r')] 

+  i[X(r')£j.,Y(r)-Y(r  )f|./X(r')]|^.  (43) 

Expanding  out  the  product,  noting  that  the  ensemble  average  of  an  odd  number  of 

2  2 

terms  is  zero  for  gaussian  statistics,  that  X  +  Y  is  the  intensity,  and  that  X  and 

Y  are  statistically  independent  and  distributed  identically  yields 

=  5  Yp.  Yr  <I(l)I(r')>  +  2<X(r)X(r')>£p.  Yr'<Y(r)Y(r')> 

-  27p,<X(r)X(r')>  •£p<Y(r)Y(r')>.  (44) 

(9) 

Fante  has  shown  that  the  intensity  correlation  is 

<I(r)I(r')>  =  1  +  exp-D^(|  r-r'l )  +  TAIL.  (45) 

We  drop  the  taQ  term  which  is  small  (see  appendix  B).  The  correlation  of  the  X,  or 

Y  components  of  is 

<X(r)X(r')>  =  <Y(r)Y(r')>  =  \  exp[-  \  D^(|r-r'|)].  (46) 

The  ensemble  average  of  ^becomes 
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(47) 


^  =5 


which  is  identical  to  Eq.  40  and  completes  the  proof.  In  the  next  section  we 
consider  a  simple  connection  formula  between  the  weak  and  strong  turbulence 
limits. 


C.  <0^>  CONNECTION  AND  ENGINEERING  FORMULAS 


We  now  have  an  equation  for  both  the  weak  and  strong  turbulence  regimes, 
given  by  Eqs.  36  and  41.  respectively.  We  will  connect  these  formulas  smoothly 
and  accurately  to  order  (i.e.,  to  order  <Tj)  as  follows 


<el>  . 


11/3 


1+z 


11/3 


400  c^^^z 


x2/^M^(x)|l+z^“'^(l-(£x)^/^zJexp(-2(€x)^/^z^]|dx 

2  ’ 

9x(kD)^(l-«^)  (1+z  (48) 

n 


where  A  and  B  are  defined  in  Eqs.  36  and  41,  respectively.  We  have  used  this  2- 

axis  equation  to  calcidate  the  focal  plane  1-axis  root  mean  square  (rms)  wander 

angle  of  a  transmitted  beam.  The  results  are  shown  in  Figures  4-33.  The  curves 

show  the  rms  wander  angle  in  microradians  vs  the  propagation  range.  Each  figure 

is  for  a  particular  wavelength,  strength  of  turbulence,  and  obscuration  ratio,  as 

indicated  on  the  figures.  Each  figure  contains  four  curves  each  of  which  is  for  a 

different  aperture  diameter.  Figures  4-27  are  for  diameters  of  0.1,  0.2,  0.3,  and 

0.4  m  while  Figs.  28-33  are  for  0.25,  0.5,  1.0,  and  1.5  m  with  the  smidler  diameter 

curves  lying  higher  in  the  figures.  Each  curve  is  qualitatively  the  same  with  the 

1/2 

rms  angle  increasing  as  z^  ,  reaching  a  maximum  then  falling  to  a  constant  value 
depending  (weakly)  on  the  wavelength,  obscuration  ratio,  and  aperture  diameter. 
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An  engineering  fcrmula  suitable  for  hand  calculation  has  been  determined 
from  Eq.  48  which  is  valid  for  obscuration  ratios  up  to  approximately  0.5.  Using 
the  asymptotic  limits  given  by  Eqs.  36  and  42,  gives  an  engineering  formula  for  the 
total  2-axis  mean  square  angle  of  arrival  as 


where 


<ol>  = 


400€®^^z  A(d) 
n 


(1-5^)  9ir(kD)^[l+z 


1  + 


,.20  5/3  .... 

z^A(<J) 


Aid)  =  [0.27534(1+ -  0.6 x  ^ 


o  +  4x^sin^a 


X  [cot(a-9)  -cotolldx 


S  0.275  -  0.989 


2 


(49) 
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RMS  Wander  Angle  of  a  Transmitted  Laser  Beam  vs  Rang 


RMS  Wander  Angle  of  a  Transmitted  Laser  Beam  vs  Range 
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Fig.  12.  RMS  Wander  Angle  of  a  Transmitted  Laser  Beam  vs  Range 
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Fig.  17.  RMS  Wander  Angle  of  a  Transmitted  Laser  Beam  vs  Range 
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Fig.  26.  RMS  Wander  Angle  of  a  Transmitted  Laser  Beam  vs 
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Fig.  29.  RMS  Wander  Angle  of  a  Transmitted  Laser  Beam  vs  Range 
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V.  ANGLE  OF  ARRIVAL  OF  A  POINT  SOURCE 


In  order  to  determine  the  angle  of  arrival  of  a  spherical  wave  in  a  receiving 
aperture  due  to  a  point  source  in  the  laser  target  plane  it  is  necessary  to  calculate 
the  ensemble  average  as  given  in  Eq.  25.  Making  use  of  Eq.  27,  we  note  that  the 
required  ensemble  average  is 


r(r  r  r  r  )  =  /  \ 

nil,  12, 13,  r^;^  \  (At  /.ut=:v.t.2„,2  A 


(/Uo(£)u;(£)d2p)2 


(50) 


Consider  the  integral  of  the  intensity  over  the  aperture.  For  weak  turbulence  the 
intensity  does  not  fluctuate  much  from  the  vacuum  limit  and  the  integral  can  be 
taken  outside  the  ensemble  average.  As  the  strength  of  turbulence  increases  the 
amplitude  coherence  length  decreases  with  its  maximum  value  being  approximately 
equal  to  ^z^/k  and  decreasing  to  a  value  of  the  order  of  the  phase  coherence 
length  p  (z).  For  the  range  of  ap«*ture  sizes  we  are  interested  in,  there  are  many 
amplitude  correlation  patches  within  the  aperture.  The  central  limit  theorem 
implies  that  the  integral  of  the  received  intensity  is  then  a  normally  distributed 
random  variable.  By  noting  that  is  normally  distributed  under  strong  turbulence 
conditions,  an  approximate  result  may  be  found  for  Eq.  50  as  follows.  Let  x  be  the 
normally  distributed  random  variable  representing  the  integral  of  intensity  over  the 
aperture.  Then 


p(£)U'^(£)d  p 


-2 


lim 

€-«0 


//doc* 


,  -€(oc+y)  . 

d7e  'sinaxsinyx. 


(51) 


Expanding  the  sin  functions  into  exponentials  yields 


sinax  sinyx  = 


^i  ( a  +7)x+ i(  a  ■h^)x_gi  ( a  -^)x_^-  i(  a  -y  )x 


We  next  substitute  this  expansion  into  Eq.  50  and  consider  each  term  in  the 
ensemble  average  separately.  Let  P  =  +(a+y),  then  for  each  term 


lim 

WliW2,W3,W4-*0 


Since  all  the  factors  in  the  exponential  are  normally  distributed  with  zero  mean, 
the  ensemble  average  on  the  right  hand  side  of  Eq.  (52)  becomes 

exp  -i/3  <x>+  |<[wj Uj+W2U2+W3U3+W4U  * -i^(x-0^)I 

We  have  simplified  the  notation  letting  the  position  coordinate  be  indicated  by  the 
subscript  on  the  field  variable.  This  exponential  factor  written  explicitly  yields 


<  exp  I  WiUj+W2U2+W3U3+W4U4'|>exp  |  -i0<^>-l/2ph<x^>-<x>h 

X  exp  i0<(x-<x»(WjUj+W2U2+W3U3+w^l/4)>|. 

The  cross  term  between  (x-^x^)  and  the  field  variables  U  is  zero  since  the  IT's  are 
normal  and  there  are  an  odd  number  of  field  variaUes  in  the  ensemble  average, 

l.e.y<U^(£>Up  (£)U^(r.)>d^£  =  0.  This  implies  that  the  ensemble  average  shown  in 
Eq.  50  becomes 


This  means  that  for  strong  turbulence,  the  angle  of  arrival  for  the  receiver  case  is 
modified  by  a  multiplicative  constant  which  is  the  normalized  ensemble  average  of 
the  square  of  the  integral  of  the  received  intensity  over  the  aperture.  That  is. 
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For  weak  turbulence  the  resiits  for  the  transmitted  beam  and  receiver  case  are 
identical. 

We  must  now  find  <l/x  >.  We  have  already  shown  that  the  ensemble 
average  of  one  of  the  exponential  factors  in  the  expansion  of  the  sin  terms  of  Eq. 
(51)  is  given  by 

exp|-i0<x>-i02j^xV<x>2]|, 
where  0  Now  let 

a  =  <x>  =  y U^(£)UJ'(£)  d%  =  (l-6\ 

and 

a^F  =  <xb-<x>^ 


=  I^expl-D^(|£r£2l» 

=  /xM^(x)exp[-2€®^^Zj^x®/^]dx.  (55) 


For  very  strong  turbulence  and  d  ^0.7  an  asymptotic  expression  for  F  can  be  found 
by  approximating  M^(x)  =  IVI^(O)  =  1-d^.  Then 


F 


^  48r(6/5) 
5x(l-d2) 


(56) 


and  upon  taking  the  ensemble  average  of  Eq.  51  we  obtain  that 


00  00 

<l/x^>  =  Wm  f  J  dadTre~*^“^^ 
«— 0  .00  00 


sinax  sin7X 


=  -  lim 
€-*0 


i// 


2  Q  Q  dadTr^exp[-€(o[+7)]  j  exp[-0.5(a+7)^a^Fl 

X  cos  (a+7)a  -  exp[-0.5(a-y)\^FI  cos(a-7)aj^ 

For  the  term  with  a  +  7  alone  we  change  variables  by  letting  0 
integral  becomes 


n  Vi  vx  ^ 

^  22 

=  lim  /  ofda  e  ®  ^^cos  aa 


™  r. 

p 

=  /  ada 


-a^Fa^/2 
e  cosaa 


-a^Fa^/2. 
a  e  da 


The  second  term  is  treated  in  a  similar  manner  by  letting  0  =  7-a. 
becomes 


lim  lae"^^" 

=  -lim  J  tinf‘ 


f 


cos/Ja 


lim  j(  dae  ^®e  ®  ^^cosaa +2c ae  ® 
€— 0_ 


-  /  rf  -a^Fa' 

-  -  (p  adae 


72 


cosaa 


(57) 

a  +  7  and  the 


This  term  then 
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Therefore 


<l/x2>  =  -i- 
a^F 

_  _1_ 
a^F 


00 

1-a  J"  si 


-a^Fa^/2 

sinaa  e  da 


00 


(4) 


k=0  (2I<+1)!! 


(58) 


Note  that  an  asymptotic  solution  for  can  be  found  when  F  is  very  small  by 

considering  the  integral 


lim 
6-0  0 


J"  ada  cos  a 


-£a  -a^Fa^/2 
a  e  e 


and  expanding  the  term  containing  F  as 

c^>  -  lim  ^  2n+l^ 

c  >  -  lim  ^  rrl — o")  n  “ 


<l/x" 


£  —0  n=0  ^  ^ 

00 


cosaa  e 


~€a 


_  1  'y  (2n+l)!/F\"^  1  „  .,p.,,p2.  .  .  , 

■  i^O  --3H*3F*15F  -  I 


(59) 


Finally  we  have  that  the  correction  factor  for  the  receiver  case  is  given  by  Eq.  58 
multiplied  by  a  .  We  have  shown  that  F  behaves  as  z  ^  for  z  »  1.  Since  we 
have  neglected  contributions  from  the  tail  in  the  intensity  correlation,  which  vary 
as  z^  in  the  determination  of  it  would  be  consistent  to  set  the 

correction  factor  to  unity.  What  is  interesting  is  that  the  correction  factor  is 
greater  than  unity  for  small  F  then  decreases  below  1  for  larger  values.  A  plot  of 
the  square  root  of  the  correction  factor  vs  F  is  shown  in  Fig.  34.  Note  that  the 
approximation  breaks  down  when  the  central  limit  theorem  is  not  valid,  which 
probably  occurs  f  or  F  values  greater  than  about  0.2. 
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VL  CORRELATION  BETWEEN  CENTROID  AND  TILT  ANGLE  OF  ARRIVAL. 


Besides  using  the  centroid  method  to  define  angle  of  arrival  it  is  possible  to 

define  the  angle  of  arrival  by  minimizing  the  mean  square  phase  error  across  the 

(8) 

aperture,  as  previously  done  by  Fried.  This  tQt  emgle  is  defined  for  the  x- 
direction  in  weak  turbulence  by 

/x0(r)W(r)d^r 

)^r'^W(r)dr 

where  0(r)  is  the  phase  as  measured  by  interferometric  techniques  and  W  is  the 
aperture  function  for  the  intensity  rather  than  the  field.  These  two  methods  can 
be  compared  by  calculating  the  corresponding  correlation  function,  i.e.. 


1/2 


where  we  have  already  defined  the  centroid  wander  angle  (see  Eq.  24)  as 


4 


=  -i_  I. 


W(r)Up(r)[7pUjr)W  (r)]d  r 


^W(r)|^I(r)d^r 


Assuming  a  uniformly  illuminated  aperture  we  have 


2  y 

v<%>  = - ^ 


and 


1/2 


j^j[y‘(R^V/4)D^((£i)dWRj^^^  ’ 

-iy/2J\JMVV%)6\ 

»  o  — r  0  — 


172  • 


(61) 


-67- 


T 


Therefore  the  correlation  function  is  (setting  =  1) 


-i  2  ff[r’  £.u: 

EE“ 


tr'  £rU^(r)l«(r')U^(r)  d^rd^r' 


Vp  'JJ(R^-p^/4)D^(|p|)dWRl^^^r  r^xdx  \ 


°  /.  \5/3 


•(62) 


For  weak  turbulence  condition  we  assume  that  ^0  =  0^^=  using  the 

results  of  Ref.  11  we  obtain  that  '  ^ 


-  i\l2  //r'. 


= 

P  c 


XL 


i  ip-Q'(r)0(r')U^(r)  d^rd^r' 


fxD^*5/3^J  /xWW(x) 

]l/2 

^  0 

(63) 


where  is  given  in  Eq.  15.  For  weak  turbulence  conditions,  the  ensemble  average 
of  ^(r  )UQ(r)  rpUQ(r)  can  be  shown  to  be  equal  to )  by  methods  similar 
to  that  used  in  section  IV -A.  Hence 


'9  6^  = 


il 

'i\j2  .2 

r-r))d  rd  r 

5/3irD^C®/\^[N^ 
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J 

^  x^^^dxM^  (x)dx 

1/2 

.  / 


(x)dx 


0 


1/2 


V  V’ 

For  5-0,  the  integrals  in  Eq.  64  can  be  performed  analytically  yielding 


(64) 


C.  ^  =  0.9919. 

'^p  c 

A  plot  of  Cg  ^  vs  the  obscuration  ratio  is  given  in  fig.  35.  Note  that  a  high 
degree  of  corrfta^ion  exists  for  all  values  of  5,  i.e.  at  5  =  0.7,  ^  is  still 

=  0.96  02.  P  c 
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For  strong  turbulence  conditions  z  »  Zg,  this  analysis  is  not  valid  since 
"phase"  is  not  a  measurable  quantity.  Phase  is  determined  analytically  from  a 
diffraction  or  intensity  measurement.  On  the  other  hand  the  angle  defined  by  the 
centroid  is  a  measurable  quantity.  Once  a  method  of  tat  determination  by  a  least 
squares  fit  of  phase  over  an  aperture  is  defined  in  terms  of  measurable  quantities 
we  will  be  able  to  complete  the  analysis  under  strong  turbulence  conditions. 


r 


Vn.  ENGINEERING  FORMULA  FOR  LASER  BEAM  PROFILES 

It  has  been  shown^^®^  that  an  engineering  formula  giving  approximate  results 
for  the  far  field  pattern  of  a  uniformly  illuminated  aperture  without  tilt  correction 
can  be  written  as 


1(0)  =  TI  R  exp(-0^/0^) 
o  o 


(65) 


where  T  is  the  atmospheric  attenuation  factor  due  to  absorption  and  scattering,  I 

-1  ^ 
is  the  peak  vacuum  irradiance,  0^,  the  1  sigma  (e  )  ang^ar  divergence  of  the 

beam  is  given  by 


2 

0  =  0 
''o  vac 


♦  d2/[2pJ(z)1  ♦ 
2 


(66) 


2  2 

0j  is  the  2-sigma  rms  jitter  angle,  R  =  ®nd  m  is  the  aperture  diffraction 

limited  performance  factor.  Under  weak  turbulence  conditions  it  has  been 
shown^^^^  that  the  short  term  average  (i.e.,  tilt  corrected)  beam  profile  is  similar 
to  the  beam  profile  in  the  absence  of  turbulence.  The  difference  is  that  the 
smaller  turbulence  scales  cause  wide  angle  scattering  with  the  result  that  the 
central  lobe  decreases  in  value  compared  to  its  vacuum  value.  This  decrease  in 

intensity  can  be  described  heuristically  by  a  factor  e(-^D)»  where 

2-2  5/3  ** 

0D  =  2.51  X  10  €  z  .  It  is  also  known  that  as  the  turbulence  strength  increases 

o  n 

that  the  short  term  average  intensity  approaches  the  corresponding  long  term 
average  value.  For  this  reason  we  write  a  general  formula  for  the  tilt  corrected 
far  field  or  focused  laser  beam  profile  of  a  uniformly  illuminated  aperture  as 


h’C 
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TC 


(67) 
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l+f,e  z^A(«) 


and  A(i)  is  given  in  Eq.  49.  Note  that  Eqs.  66  and  67  contain  the  same  symbol  for 
the  effect  of  laser  platform  jitter  on  the  system  performance.  These  are  not  the 
same  quantities  since  tat  cwrection  can  correct  for  a  great  deal  of  platform  jitter. 
The  quantity  6j  in  Eq.  67  should  be  considered  as  the  uncorrectable  or  residual 
jitter. 

In  this  report  we  have  been  concerned  with  the  effectiveness  of  tilt 
correction  as  the  strength  of  turbiilenee  becomes  more  pronounced.  We  have  found 
that  as  the  strength  of  turbulence  increases  the  mean  square  an^e  of  arrival  of  a 
transmitted  beam  increases  linearly  for  z  <  Zg  reaches  a  maximum  (z  ~  Zg),  and 
then  decreases  for  z  >  Zg  to  a  constant  which  is  of  the  order  of  the  diffraction 
limited  spot  size  of  the  aperture.  This  means  that  the  effectiveness  of  tilt 
correction  decreases  markedly  for  z  ^  Zg.  A  measure  of  this  effectiveness  could 
be  found  dividing  the  laser  beam  mean  square  tot  an^e  by  the  long  term 
average  mean  square  angular  divergence  of  the  laser  beam.  When  this  ratio 
becomes  very  small,  tilt  correction  is  not  very  useful.  We  have  also  shown  that  a 
phase  measurement  is  sufficient  to  determine  the  tilt  correction  under  conditions 
of  weak  turbulence;  however,  under  strong  turbulence  conditions  it  may  be 
necessary  to  measure  the  tilt  of  a  received  spherical  wave  directly  in  order  to 
perform  the  reciprocal  correction.  This  direct  measurement  would  require 
measurement  of  the  location  of  the  centroid  of  intensity  in  the  focal  plane  of  the 
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receiving  aperture. 

In  this  report  we  have  given  simple  engineering  form liae  for  determination 

of  mean  square  angle  of  arrival  (Eq.  49)  and  the  relative  peak  i''‘ensity  of  a  tit 

corrected  laser  beam  (Eq.  67).  For  example,  consider  the  oa^-j  of  fairly  strong 
9  _1 4  _2  /3 

turbulence  vnth  C  =  10  m  and  X  =  0.5  gm.  The  saturation  range  is 

approximately  4  km,  and  p  (zj  =  0.75  cm.  Using  an  unotacured  aperture  diameter 

of  0.25  m  or  (4  =  33)  and  a  propagation  range  of  10  km  the  2-<t  root  mean  square 

wander  angle  is  3.48  prad  and  I^j  =  5.8  x  10"^,  i.en  no  tilt  correction  is  possible  for 

this  case.  At  a  range  of  2  km  the  2-aroot  mean  square  wander  angle  is  7.94  prad 
-3 

and  =  6.1  x  10  which  represents  about  a  30%  improvement  over  the 
uncorrected  tilt  case. 

This  improvement  factor  does  not  seem  sufficient  based  on  our  previous 
resiits  in  the  weak  turbulence  regime.^^^^  In  deriving  the  long  term  average  beam 
spread  we  assume  that  the  peak  on-axis  intensity  is  directly  related  to  the  beam 
spread,  i.e.,  that  a  gaussian  proffle  resiits  under  turbulent  conditions.  From  this 
we  have  found  that  the  turbulence  induced  beam  spread  additive  factor  is 

C  '  »v.o 

This  engineering  formiia  has  proved  very  effective  in  the  past,  for  small  values  of 
D/p^(z)  (<10),  in  determim'ng  long  term  average  spread  for  general  aperture 
distributions  which  were  gaussian  and  truncated.  If  we  use  the  results  for  the  peak 
on-axis  intensity  for  a  uniformly  illuminated  aperture  for  large  values  of  e,  we 
obtain  the  asymptotic  value 

lpa4.4072[,2fe)/D=']lo 

where  I  is  the  peak  irradiance  without  turbulence.  This  means  that  we  coiild 

®  A  «  ye 

modify  Eq.  67  such  that  ^  z”'  would  be  divided  2.204.  for  the  second  case 
g^ven  above,  =  0.5, 1^^^  would  then  be  approximately  =  0.03  or  an  improvement  by 
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a  factor  of  75  over  the  uncorrected  tilt  case. 

This  process  of  subtracting  the  mean  square  wander  an^e  from  the  long 
term  average  beam  spread  deserves  further  examination.  We  have  shown  pre- 
viously^^®^  fcr  weak  turbulence  conditions,  that  if  the  instantaneous  phase 
aberrations  of  the  atmosphere  are  cwrected  to  first  order,  i.e.,  for  tilt,  defined  by 
a  least  fit  of  the  measured  phase  of  the  return  s{*erical  wave  over  the 
transmitting  aperture,  then  the  transmitted  beam  profile  looks  like  a  vacuum  beam 
profile  with  a  gaussian  halo  caused  by  wide  an^e  scattering.  The  intensity  profile 
of  the  long  term  averaged  beam  is  gaussian  however;  therefore,  the  two  beams, 
corrected  and  uncorrected,  do  not  look  the  same  €ind  the  wander  an^e  should  not 
be  subtracted  from  long  term  beam  spread  to  get  the  appropriate  short  term  beam 
radius  or  peak  on-axis  irradiance  for  the  tilt  corrected  case.  In  order  to  obtain  the 
engineering  formula  67  we  used  the  engineering  equation  obtained  from  examining 
the  tilt  corrected  laser  beam  profile  fcr  weak  turbulence  (Ref.  10)  and  extrapo¬ 
lating  the  results  into  the  saturation  regime  saying  that  eventually  the  results 
should  retirn  to  the  long  term  average  results  as  z  »  z^.  This  process  is  extremdy 
suspect  and  must  be  examined  in  a  more  critical  manner. 
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APPENDIX  A 

VERIFICATION  OF  THE  NEGLECT  OF  THE  (BJ)^  AND 

/  9  ^ 

(SX)^  TERMS  IN  WEAK  TURBULENCE 

^2  ^2 

In  this  appendix  we  will  verify  the  neglect  of  the  terms  (Bw)  and  (SX)  in 

2  ^ 
comparison  to  7  D  in  Eq.  30  under  conditions  of  weak  turbulence.  To  this  end  we 

2 

note  that  Tatarskii  has  defined  structure  function  terms 


D,(o)  = 


-5/3 


X 

jFj  (-5/6,1,-  J-1 


A-1 


12 


where  ^F^  is  the  confluent  hypergeometric  function  . 
andXS  structure  functions  are  defined  in  terms  of  Dj,  D2 


D^(p)  =  Dj(p) 


Dy(p)  =  \  (Dj(p)-Re  D2(p)] 


=  \  ImD2(p) 


and  we  define  the  parameters 
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The  wave,  log-amplitude, 
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A-3 
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2 

where  X  is  the  inner  scale.  Noting  that  the  structure  function  is  related  to  the 

correlation  function  by  D(p)  =  2(B(0)-B(p))  and  that  the  results  for  the  structure 

functions  for  a  spherical  wave  are  related  to  that  for  the  plane  wave  (Eqs.  A-1  and 

A>2)  by  replacing  p  by  sp  and  integrating  on  s  from  0  to  1  we  can  proceed  by 

considering  different  regimes  on  p.  Fimt  we  are  generally  interested  in  the  case 
2 

when  XL  x>X  .  Since  the  inner  scale  is  of  the  order  of  1  mm  and  for  X  between 
o 

0.3  and  10  pm  and  10  cm  <  L  <10  m  this  condition  is  satisfied.  For  g  «  1  or  p 
smaller  than  the  inner  scale 
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Now  for  this  regime 


Therefore 


Similarly 


=  1.09 


r  T  \  0  071P  2.  13/6j^5/6 

[bxsw]^  =  '  " 


81B^(p)]^ 


^  0.126  <T* 


'0 


/2.x2V/« 

br) 


0 


;2irx2\'^/® 
9xl0^a2  — o 
‘'TV  XL 


«  1 


A-8 


A -9 


A-10 


A-11 


A-12 


-76- 


^2  ^2 

For  the  inner  scale  region  the  neglect  of  the  extra  terms  (^)  and 
therefore  justified. 

Consider  now  the  region  \^<  p<\/ljk.  For  p  >  X  we  know  D./,(p)  = 

i5  /3  ^  o  V' 

2[p/p  (z)r  and  we  have  shown  that 


A-13 


(2) 

For  p  between  the  inner  scale  and  the  Fresnd  radius  Tatarskii  has  shown  that 


'  2 

F  or  this  reason  )  = 
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In  order  to  determine  B^g  we  must  expand  the  confluent  hyper  geometric  function 
given  in  Eq.  A-1.  Expanding  Eq.  47,29  of  Ref.  2  we  find  for  g/Q<l  that 
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Using  Eq.  A-16 
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Finally  we  consider  the  case  of  p  >>/L/k.  For  this  case 
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where  we  have  used  Tatarskii's  definition  of  b^.  Then 
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where  [c]  has  been  defined  by  Eg.  A-16.2.  Using  the  asymptotic  expansion  for 
F  (-11/6, -11/6, -ig/Q)  as 


F  (-  11/6,  -  11/6,  -  ig/Q)  =  1  +  i(ll/6)2(Q/g) 

.  (ll/6)^(5/6)^(l/6)  +...  ^_22 

the  asymptotic  value  of  D^g  is 

Pxs  ^  -|[clC^k^(ll/6)^(5/6)^(^)^(XL)V^^ +coretant  A-23 
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which  completes  the  analysis  and  verifies  the  neglect  of  (b’ 
from  Eq.  30  for  all  regions  of  p.  ^ 


i)  terms 


APPENDIX  B 

THE  EFFECT  OF  THE  FOURTH  ORDER  CORRELATION  TAILS  ON 

Recall  that  the  calculation  of  the  mean  square  wander  an^e  involved  the 
determination  of  the  ensemble  average  of  the  intensity  correlation  function 
<I(r)I(r')>.  Fante  has  derived  an  iterative  solution  to  this  correlation  function 
under  strong  turbulence  conditions.  To  the  second  order  in  the  iteration  procedure 
the  intensity  correlation  is 


<I(rj)I(r2)> 

where 


=  1 


+  exp  -  D^Cjrj-rgl 


(^2)2/5 
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g(W) 


r  f° 

=  0.27  J  dt  r®''^l-cos  t)  J  dse'^J  (2.43  t 


8/3s-3/ll^) 
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=  1.19  qI  /  r^^/®sin^, 

X  exp  |-djt®''V'^®[4.26-2.66y]|  dt 


B-3 


Now  R  =  0.226  ex,  i.e.  g(W)  is  independent  of  z  .  Therefore,  the 

contribution  of  g(W)  to  must  go  like  or  z^"  ^^'^^^where  (<rp~  is 

the  preceding  factor  ot  g  in  Eq.  B-1.  For  large  z  this  contribution  can  therefore 

n 

be  ne^ected.  The  contribution  due  to  Eq.  B-3  is  much  harder  to  determine.  First 
take  the  gradients  shown  in  Eq.  32. 
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Note  that  x  is  the  normalized  aperture  integration  variatde.  Change  variables  by 
setting  w  =  then 


T  = 


2^  7/15  ,  00 

^  -  -1/3^  /*  -4/5  .  .  2| 

y  dy  #  w  dw  SI  n 
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For  the  uniformly  illuminated  aperture,  the  contribution  from  T  is  found  by 
performing  the  double  integration  over  the  aperture  whidi  yields 


^  IP  * 


7/15  r  -1/3  .  /^-4/5  ^  .2 

=  1 2. 63  €  z  /  y  dy  /  w  dw  si  n - 
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Divide  the  w  integration  into  two  regions  which  are  less  than  or  greater  than 

c  yO  1  1  /I  R 

(5.52/<)  '  .In  the  latter  region,  the  Bessel  function  term  is  less  than 

unity  and  the  exponential  factor  which  is  smaller  than  exp  -  const  7}^^^  where 

const  =  27.59  t  The  factor  e  const  goes  to  zero  very 

rapidly  with  and  can  therefore  be  neglected. 

The  remaining  integral  is  divided  into  three  terms  which  depend  on  the 

range  of  y  integration.  These  are  0<y 

(w/(72)^/^<y  <€w^^^/5.52  and  €w^^^/5.52z^^''^®<y  <  1.  For  all  ranges 

the  exponential  is  replaced  by  exp(-1.6  w).  For  the  first  range  the  sin  function  can 

12/5  2  2  12/5 

be  replaced  by  1/2  while  in  the  last  two  ranges  it  is  replaced  by  w  /y  . 

In  the  first  and  last  range  the  Bessel  function  factor  is  replaced  by  unity  while  in 

the  middle  range  it  is  replaced  by  the  asymptotic  approximation  for  large 
3 

argument,  I./ttQ  .  For  the  first  range 
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This  term  decreases  like  and  can  therefore  be  neglected.  The  contribution 
from  the  third  range  is 

_  0.453  /  8/5  -l/6w  f  -7/3. 

0-'  ^  ^ 


2  yP* 
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dll  =  €w2/5/(5.52z^^^^^®) 
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This  term  goes  like  z^"  '  .  Therefore  it  can  be  dropped.  The  middle  range  gives 
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i.e.,  this  term  also  goes  like  z^^”  '  and  can  also  be  dropped.  This  concludes  our 
proof  that  the  tail  terms  can  be  neglected  in  the  determination  of  the  mean  square 
wander  angle  of  a  transmitted  laser  beam. 


REFERENCES 


1.  The  whole  March  1977  edition  of  J.  Opt.  Soc.  Ann.  is  devoted  to  various 
adaptive  optics  techniques,  both  experimental  and  theoretical. 

2.  V.  1.  Tatarskii,  "The  Effects  of  the  Turbulent  Atmosphere  on  Wave  Propa¬ 
gation",  National  Technical  Information  Service,  Springfield,  Va.,  1971. 

3.  R.  F.  Lutomirski,  R.  E.  Huscke,  W.  C.  Meecham  and  H.  T.  Yura,  "Degrada¬ 
tion  of  Laser  Systems  by  Atmospheric  Turbulence",  Rand  Corp.  report 
number  R-1711-ARPA/RC,  June,  1973. 

4.  Capt.  J.  Taylor,  AFWL/LEAPS,  private  communication. 

5.  R.  E.  Hufnagel,  "Variations  of  Atmospheric  Turbulence",  Paper  WAI  appear¬ 
ing  in  the  Digest  of  Technical  Papers,  Topical  Meeting  on  Optical  Propaga¬ 
tion  Through  Turbulence,  ;July  9-11,  1974,  Univ.  Ccdorado,  Boulder,  Colo. 

6.  M.  Tavis,  Private  Communication. 

7.  S.  F.  Clifford,  G.  R.  Ochs,  and  R.  S.  Lawrence,  "Saturation  of  Opticed 
Scintillation  by  Strong  Turbulence",  jour.  Opt.  Soc.  Amer.,  M,  148,  (1974). 

8.  D.  L.  Fried,  J.  Opt.  Soc.  Am.,  56,  1372  (1966). 

9.  R.  L.  Fante,  "Irradiance  Scintillations:  Comparison  of  Theory  with  Experi¬ 
ment",  AFCRL-TR-74-0626,  20  Dec.  1974. 

10.  H.  T.  Yura,  "Atmospheric  Turbulence  induced  Laser  Beam  Spread",  Apl. 
Opt.^,  2771  (1971). 

11.  M.  T.  Tavis  and  H.  T.  Yura,  "Short  Term  Average  Irradiance  Profile  of  An 
Optical  Beam  In  A  Turbulent  Medium",  Apl.  Opt.  1^.,  2922,  (1976). 

12.  M.  Abramowitz  and  I.  A.  Stegun,  "Handbook  of  Mathematical  Functions, 
N.B.S.  No.  55,  Dec.  1965. 


-85- 


LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory;  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution:  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon:  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields:  space  astronomy,  x-ray  astronomy:  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  apace  systems. 
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